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Abstract 

The goal of this paper is to present a novel approach for total varia- 
tion regularization and Sobolev minimization, which are prominent tools 
for variational imaging. Thereby we use derivative free characterizations 
of the total variation semi-norm and Sobolev semi-norms of functions re- 
cently derived by Bourgain, Brezis, Mironescu and Davila. Their analysis 
is to approximate the semi-norms of a function by singular integral opera- 
tors. With this characterization we derive a series of novel regularization 
methods for total variation minimization which have as a novel feature a 
non-local double integral regularization term. 



1 Introduction 

Given noisy image data f s , total variation denoising (see |14p consists in min- 
imization of the functional 

^(/) :=\jtt-f) 2 {x)dx + a\Df\. 

The minimizer is a smoothed approximation of / . In the above functional \Df\ 
denotes the total variation seminorm of / and a > is a positive constant. The 
first summand of the functional above is called the fidelity term and penalizes 
the deviation of an image / from the data f s . The second summand is named 
regularization term and penalizes the rate of change within /. While the first 
summand provides that the outcome of the minimization process (the denoised 
image) preserves similarity to f s , the second term is intended to reduce the 
oscillations within the argument / in order to generate an approximation to f s 
that is free of the inherent noise. In fact, this strategy has proven to be successful 
and even more, the total variation seminorm has proven to be superior to other 
regularization terms penalizing the rate of change in the sense that edges within 
the image are preserved and not blurred. 

Another well-known choice for the regularization term is for 1 < p < oo the 
p-th power of the Sobolev (l,p)-seminorm | • |i. p given by 

\f\i, P ■■= (7 l V /(*)| P dx 
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Exchanging it with the TV-scminorm in the above functional yields the func- 
tional 

F p (f) -=\j (f-f) 2 (x)dx + a\f\l p . 

In this paper we are concerned with a analytical approach to approximate 
the functionals T v in a variational sense. Recently, new derivative free char- 
acterizations of the Sobolev spaces W 1,p and the space BV of functions of 
bounded total variation function have been obtained by Bourgain, Brezis and 
Mironescu [3] and Davila [9] - This work has been refined and supplemented 
by Ponce in [13j . These characterizations provided new derivative-free approx- 
imations of the (l,p)- and total variation seminorms that are obtained by ap- 
proximating the respective seminorm of a function / by double integrals over 
the p-th power of the difference quotient function of / multiplied with weighting 
kernel functions that form an approximation the identity. 

In detail, let for 1 < p < oo and any measurable function / 

KAf)--= I I ~ f ^ )l \ n (x-y)dxdy. (1) 

JnJn \ x ~ V\ 

The functions tp n are non-negative, radially symmetric, and radially decreasing 
functions from L 1 (M Ar ) satisfying that for every 5 > 



n — >oo 



lim / tp n (x) dx = (2) 

"" l{x:\x\>S} 



and for all n E N 

tp n {x) dx — 1 . (3) 



Conditions ([2]) and ([3]) imply that the unit mass of the functions ip n concentrates 
around the origin as n strives to infinity. 

If 1 < p < oo, then by [21 Theorem 2] there exist real constants K p ^m such 
that for every measurable function / defined on Q, 

lim Tllif) = [ \Vf(x)\ p dx = \ff , 
K-p,n Jn 

where the limit is +oo if / does not belong to W 1:P (il). For the space of 
functions with finite total variation a similar characterization with p = 1 in the 
integral expression holds true [9]: 

! lim n l n (f) = \Df\ . 



These approximations of the seminorms give rise to the following approxi- 
mations of the functionals J- p from above. Let for 1 < p < oo 

■= \ I (f - f) 2 (x) ^ + ~^-ni{f). (4) 

In this paper we give a variational analysis of the approximation of the 
functionals T v by the functionals for n — > oo. In particular, we will show 
that all of the approximating functionals J 7 ? have a unique minimizer /„ and 
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that the sequence of minimizers (/„) has an accumulation point / that is indeed 
the unique minimizer of the limit functional T v . That is, the minimizers of J 7 ? 
approximate the minimizer of T v . Most of these results are applications of 
results found in [3], [2] and [TB"] . 

Further, we show how the regularization functional 7£„ evaluates for some 
important examples. In this way, taking into account the above mentioned re- 
sult of the paper, we derive a series of numerical schemes for total variation 
minimization. In fact well-known numerical methods for total variation mini- 
mization can be derived, as well as a series of new ones. In future work this 
approach could perhaps be used to clarify the relationships between existing 
and seemingly unrelated numerical and analytical approaches or to supply an 
analytical basis for existing numerical schemes. A particular consequence of our 
approach is that by the used approximation it turns out that total variation 
minimization can be considered a bilateral filtering 16J. There is still an ongo- 
ing discussion on comparing qualities of total variation denoising and bilateral 
filtering numerically (see e.g. [H])- The results of this paper shed some addi- 
tional light on this topic from an analytical point of view. Another issue is that, 
in contrast to total variation minimization, the new functionals do not require 
the differentiability of the total variation measure, and thus the derivation of 
numerical schemes can be considered in a function setting. 

Notations. We summarize some further general assumptions and notations 
that will be used throughout the rest of this paper. The set C M. N is a bounded 
open set with C 1 -boundary. The symbol Tt N ~ 1 denotes the (N— l)-dimensional 
Hausdorff- measure in M. N . The symbol C denotes the A^-dimensional Lebesgue 
measure in K . 

For a real function / as usual /+ := max(/, 0) and /_ := — min(/, 0) sym- 
bolize the positive and negative parts of /, respectively. 

The letter p is always used as integration index and always satisfies 1 < p < 
oo. Sometimes the range of p will be further restricted. The letter p* denotes 
the dual index to p and is given by 

i + l-i 

p p* 

if p > 1 while p* = oo if p = 1. We denote by ||/|| p the L p -norm of / on il. The 
space of i p -functions on f2 with mean value zero is symbolized by 



The characteristic function of a set A C R is denoted by \A- 

Let C£°(fi) be the space of infinitely differentiable functions from f2 to R 

with compact support. For the Tl^-seminorm of a locally integrable function / 

we use the symbol 




The mean value of an integrable function / on Q is symbolized by 




\Df\ = sup 
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The space of functions of bounded variation on fl is the set 
BV(Sl) ^{/S^fl) : \Df\ <oo}. 
The constants K p ^ are defined by 

T^-i^-i) ! s h-x \{e,a)\ P dH N ~\a) if N > 1, 
1 i£ AT = 1. 
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The radial and radial decreasing functions </?„ give rise to the monotone 
decreasing functions (p n on (0, oo) defined for all r > by 

lfin{r) := <fin{x) 

where x is a point in R w with |a;| = r. 
Let 

S:L l {VL) -> RU{+oo} 

/ - ^II/-/ A ||2 

denote the first summand of the functional J 7 ^. The functional attains the value 
+oo if / - f not in L 2 . 



2 Variational Analysis 

Let 1 < p < oo. In the following we prove existence and uniqueness of a 
minimizcr of the functional defined in (U]). 
We begin with the following Lemma. 

Lemma 1. Let f be a real measurable function that does not belong to L p (fl). 
Then 

K(f) = 

Proof. By the hypothesis 

f \f(x)\ p dx = +oo. 
Jo, 

By the assumptions on tp n there must exist a radius r > with (p{r) > 
and ip(s) > <p{r) for all s < r. 

Since f2 is bounded Cl can be covered by finitely many balls of radius r, and 
there must exist a measurable subset M of fl with positive measure satisfying 

/ \f(x)\ p dx = +oo and diam(M) < r. 
Jm 

It follows that for any (3 € R the function / — j3 ^ L P (M); in particular 
f - f(x) LP(M) for almost all x & M. 
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Now we have 

K(f)>[ I lfi f~ f{ * )l \ n(x-y)dyd X 

Jmjm \ x — y\ 

>4 / / \f(v) - f(x)\ P Vn(x - V) dydx 

rF JM JM 

and by monotonicity 

-J f \f(v) - f{x)\*<p n (x - y) dydx >^ / / \f(y)-f(x)\*dydx. 
' Jmjm r jmjm 

But as inferred above for almost all x € M 

I \f(y)-f(x)\ p dy = +oc. 

JM 

Thus, 

/ / \f(y) - f(x)\ P dydx = +oo, 

and by the above estimations the proof is complete. □ 

We continue with the following lemma about the weak lower semicontinuity 
of the regularization functional. 

Lemma 2. Let 1 < p,q < oo. For all neFj, lZ v n is weakly lower semicontinu- 
ous on L q (Q), that is, 

ni(g) <liminf1lZ(g k ) 

k — >QO 

for every sequence (<j%) S L q (£l) that converges weakly with respect to the L q - 
topology to a function g £ L q (fl). 

Proof. By a standard result of convex analysis (see e.g. [7]) it suffices to show 
that the functional lZ p n is convex and lower semicontinuous on L q (il). Both 
properties of JZ p n are established below by representing the functional as the 
pointwise supremum of convex and lower semicontinuous functionals on L q (Q). 

Let g £ L q {yi). Note that TZ^ (g) may be viewed as the p-th power of the 
p-norm on f2 x Q of the following measurable function 

\y-x\ 

Let A := {(x, x) : x £ ft} and 

U E = {(x,y) e!lx!l: dist ((x, y), A) > e} . 

By monotone convergence 

ll$llz-(nxn) =™p{ll$ lnJ| LP( n E )} • 

Let £ be sufficiently small such that Il e is not the empty set. We use the 
notation 

B s p , := [fj £ C c (fi x O) : r?| (nxn) \ ne = and \\ij\\ LP . {Uc) < l} 



5 



for the set of continuous functions on fl x fi with compact support and LP -norm 
less than one that vanish in the complement of II e in il x ft. 

By duality, the p-norm of g on II e can be represented as the following supre- 
mum: 

< 1 



P lnJI L p(n e ) = sn P \ J n 9(x,y)v(x,y)d(y,x) : ?7eC c (n £ )and |M| L p„ (rie ) 

Using that the functions from B z Pf vanish in (ftxft) \ II £ this supremum can be 
rewritten as 



sup / / g(x,y)fj(x,y)dydx. 
■eBt. Jn Jn 



j?eBJ» Jn Jn 

We need the following property of the functions i) from the set B| t : for all 
such f] the function 

x~ f 4(x-y)j^ldy (6) 
Jn \x-y\ 

belongs to L q {Q). To see this note that for p > 1 we get by an application of 
Holder's inequality and by using the vanishing property of i) and the hypothesis 
that ip n has integral one 

/ ( I 4{x-y)f^dyX dx<[([ vkx-y)M^dy\ q dx 

Jn \Jn \y-x\ J Jn \Jn\B( x ,e) \V ~ x \ J 

<f(f <pUx-v)^dy) dx 

Jn \Jn\B( x ,e) £ ) 

<{^ym^dx 

< CO. 

For p = 1 an analogous computation can be carried out similarly. 
In the same way it is proven that for all fj £ B^ the function 

I/h / yi{x-y)f^dx (7) 

Jn \ x - y\ 

lies in L q (£l). By Fubini's theorem we can infer that 

/ / \g(yM(x-y) l -^^-dydx = [ \g(y)\ [ <p» (x - y) dx dy 

JnJn \ x -y\ Jn Jn \ x ~ V\ 

< oo 

and thus the integral 

/ [ g {y) V l{ X -y)f^Adydx 

JnJn \ x V\ 

converges. 
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Using the preceding results we can rewrite 



IISlMnxQ) =sup||3 | n J| LP(n 

e>0 v 



^sup sup / / g(x,y)fj{x,y)dydx 

iff — Tjixy y ) 

= sup sup / / g(y)ipn(x-y)-, — '—rdydx- 
E>o^ess, VJn Jo F~2/| 

50) / Vn{x-y)-, '—rdydx 

Jn \x — y\ 

Exchanging the two variables in the first integral expression followed by an 
application of Fubini's theorem shows that the difference from above is equal to 

sup sup / g(x) / ipn{x~y)- - — - dydx- 

e>o fjesj, \Jn Jn F ~ V\ 

g(x) / '—rdydx 

Jn F - 2/1 

= sup sup / / ifn{x-y) ^ 1— ! — dycfcc . 

e>0i)GBJ, Jn Jo 12/ — ^ I 

For every e > and any fj S B|^ the inner integral expression in the double 
integral expression above gives rise to the function 

\t ^l(y,x) - fj{x,y) 

fn [x - y) : : dy. 

\y ~ x\ 

Just as in ([6]) and |J7]) it is shown that all the functions belong to L q (fl). 
Thus, for all e > the functional 



g{x)hf l {x) dx 

is continuous on L q (il) for all fj £ B^, and since for all g € L q (fl) 



0^n) p (9) = llffllzxnxn) = SU P SU P / g{x)hf,{x)dx 

e>o»7eB^ Jn 

1 

the functional (T^) 13 is the pointwise supremum of continuous functionals on 
L q (Vl). We conclude that the functional IZ^ is indeed the supremum of continu- 
ous functionals on L q (fl) and hence, is lower semicontinuous on this space. The 
convexity of IZ^ can be proven in a likewise manner. □ 

Below we will use some compactness results of Bourgain, Brezis, and Mironescu 
[3] which are as follows. 

Theorem 3. Let 1 < p < 00. Assume that (g n ) is a sequence of functions 
in Le(Q) such that TZP t (g n ) is uniformly bounded. Then the sequence (g n ) is 
relatively compact in L©(f2) and has a subsequence {g nk ) converging (in the L p - 
norm) to a limit function g that lies in W 1,P (Q) if p > 1 and in BV(Q) if 
p=l. 
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In the following we apply the above results to prove existence and uniqueness 
of a minimizcr of the functional T v . 

Here wc make use of a scale space property of variational denoising al- 
gorithms, that they are grey level invariant (see 1 5J ) . We note that for all 
1 < p < oo a function /o is a minimum of T v if and only if /o — f s minimizes 
the modification of T v where f s is replaced by the function f s — J n f s . Since 
the latter function has mean value zero we restrict our attention to the case 
that the mean of f s is zero, and consequently, also the mean of the minimizer 
of T v is zero. 

Proposition 4. Let 1 < p < oo and a > and assume that f E L^(£l). 

1. Then the functional J 7 ? attains a unique minimizer f n over L^(£l) that 
also belongs to L%(fl). 

2. The function f n is also a minimizer of over 

3. The sequence of numbers (7^(/n)) is uniformly bounded over n£N, and 
the sequence f n has a convergent subsequence whose limit f is an element 
of W 1 ' P (Q) ifp > 1 and of BV(Q) if p = 1. 

Proof. We begin with the proof of the first assertion and show first that the func- 
tional J 7 ? attains a unique minimizer /„ over L\(£t). Let (g k ) be a minimizing 
sequence of functions in L^(Vl) for the functional J 7 ? such that 

lim Tl{g k ) = inf J*{g) < J*(Q) = ~ \ < oo . 

Since 

\\g k - f\\l<2Tl{g k ) 

for all k G N this implies that (g k ) is uniformly bounded in L%(CL) and thus 
has a weakly convergent subsequence in L|(f2). Let us denote this subsequence 
again by (gk) and its weak limit by /„ 6 L|(f2). 

We notice that the first summand S of J 7 ? is convex and continuous on 
L%(Q) and thus, weakly lower semicontinuous on this space [7J. From Lemma [2] 
it follows that IZ^ and thus, also J 7 ? are weakly lower semicontinuous on L|(0), 
too. 

It follows that 

mfn) < liminf Ti(g k ) = inf J*(g). 

Thus, f n is a minimizer for J 7 ?. 

As the first summand S of J 7 ? is the composition of a convex and a strictly 
increasing function it is strictly convex. By (the proof of ) Lemma [5] the second 
summand IZ^ of J 7 ? is convex and thus, we can infer strict convexity for the 
whole functional J 7 ?. The latter implies the uniqueness of the minimizcr. 

By Lemma Q] the minimizer /„ also belongs to £%.(&). 

We now show the second assertion, namely, that f n minimizes J 7 ? over L 1 (f2). 



8 



Since the mean value of f s is zero, for each g 6 L 2 (f2) we have 



S[g- g{y)dy\ = / g(x) - g(y) dy - f >(x) dx 



n 



2 



(g(x)-f s (x)) 2 dx-^g(x) dx^j 



2 



< S{g) 

and 

ni (g- J^g(y) dy^j = Tll(g). 
Thus, for every g S L 2 (il) there exists a function g £ L 2 (il) with 

m~g) < ma), 

and it follows that the function f n is a minimizer of over L 2 (Q). It is also 
a minimizer of J- p over L 1 (f2) since the first summand of equals infinity for 

jel 1 (l])\L 2 (ll). 

For the proof of the third assertion note that for all n E N 



niUn) < -mfn) < -mo) = ^-\\f s \\l 

a a la 



Therefore, the sequence {'R- P l {fn)) is uniformly bounded. By theorem [3] the 
sequence /„ has a convergent subsequence whose limit / lies in W 1,p (fl) if p > 1 
and in BV(Q) if p = 1. □ 

Notation: We denote the subsequence occurring in the third assertion of 
the previous proposition again by (/„) and use that notation for the rest of our 
paper. We denote by / the limit of (/„) for the rest of the paper as well. 

The following remark provides another justification that the mean of the 
minimizers /„ is zero. 

Remark 5. Let h be the function 

h : n -> M 
x i— ► 1. 

Note that for all g € L 2 ($1) the functional !fP is Gateaux-differentiable in direc- 
tion h and that the Gateaux-derivative (!FP)'(g;h) satisfies 

{Fl)'{g\h)= [ (g-f s )(x)dx 
Jn 

since 11? (g + Sh) = K p n (g) for all g G £ 2 (f2) and all real 6^0. Thus, we can 
infer 

"(fn-f)(x) dx = (J*)'(f n ;h) = 

2 

which means that the mean value of the minimizer /„ is equal to the mean value 
of f s that was assumed to be zero above. ■ 
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It remains to clarify whether the limit function / is a minimum of the re- 
spective limit functional T v . The concerning questions are answered to a large 
extent by a result of A. Ponce [13] in terms of T-convergence. 

We recall the definition of T-convergence in L 1 (fi) [5]. Let (F n ) denote a 
sequence of functionals mapping functions from to the set of extended 

real numbers R, and let F be a functional of this kind, too. Then the sequence 
(F n ) T-converges to F with respect to the L 1 (f2)-topology if and only if the 
following two conditions are satisfied 

• for every g £ L 1 (17) and for every sequence (g n ) in L 1 (J7) converging to g 
in the Li-norm we have 

F{g) < liminf F n (g n ); 

n — >oo 

• for every g £ there exists a sequence (g n ) in converging to g 
in the Li-norm with 

F(g) = lim F n (g n ). 



In this case we write 



We denote by sc L1 ^F the lower semicontinuous envelope of the functional 
F with respect to the strong L -topology, that is, sc^i^F is the greatest lower 
semicontinuous functional less than or equal to F. 

Ponce's result is established in a far more general setting than the one we 
are treating here. He investigates double integrals of the kind 

f\g{y) - g(x)\\ , w , 

uj : : p € [y - x) dydx 

nJn \ W — x \ J 

where u> is a continous function from [0, oo) to [0, oo) and (p e )e>o is a family 
of nonnegative functions in i 1 (R Ar ) satisfying the conditions ([2]) and ([3]). The 
functions p e are not assumed to be radial or radially decreasing. 

The functions p t induce positive Radon measures p t on the sphere S 1 ^^ 1 : 
Let B be a Borel subset of S' JV_1 and let 

R+B :={rx : r > and x £ B}. 

be the cone with its apex in the origin that is generated by B. Let 



He{B) := I p e (x) dx. 

JR + B 

The family of measures (/z e ) is bounded by 1 and thus, has a subsequence 
fi e . that converges weakly to a Radon measure p on S' Ar_1 . Moreover, let w M be 
the real function on M. N defined by 

uj^v) := / uj(\(v,(x)\) dp(a) 
Js N - 1 



for all v £ 
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Further, let u>** denote the convex lower semicontinuous envelope of the 
function uj, that is, by our assumption, the greatest convex function less than 
or equal tow. 

Finally, Ponce defines the functional 

F : -» [0,+oo] 



_ , s; ^(V ff (x)) ifgeC 1 ^), 
1 +00 otherwise. 



Ponce's result, which he proves for bounded open sets with Lipschitz bound- 
ary, is as follows: 

Theorem 6. If for all x e K w 

then 

for every g € L 1 (fi). 

Let us now apply Ponce's result to our problem. Let \i n denote the Radon 
measure induced by ip n . Since ip n is radial with integral one we get 



fJb n (B) = / ipn(x) dx = 
Jl&j-B 



n N ^ 1 (s N - 1 ) 

for all Borel subsets B of S^ -1 . 

Therefore, all the measures \i n are equal and thus, in contrast to the more 
general situation treated by Ponce, the whole sequence is converging weakly 
to a limit measure \x that itself equals all of the measures \i n . 

We note that in our case u is the function defined by 

u)(x) := \x\ p 
for all id". Thus, for all v£R N 

^(v) = J sn i \(v,a)f d»(a) = y<JV _ 1( 1 sJV _ 1) j sN i \(v^)\ P dH N ~ 1 {<j), 

and since the integrand only depends on the length of v we have 

M*) = nN^N-^ l< e ^>l P dn N -\a) = K P Mv\ p . 

where e is an arbitrary unit vector in M. N . Since both functions u> and u) v are 
convex it follows that 

(w M ) = (w )„ 

and thus, Ponce's theorem is applicable to our problem. We get 

r ~mn)-^J n f a - *) dydx = sc- L1{n) F(g). 
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for all g G ^(Q). 

We have to determine the functional sc^i^F. From theorem 1.2. in [B] we 
get that for all g e L 1 ^) 



sc z,i(n) 



laF ( $ ) = f K *M^9{x)\ P dx + f v?(j$fr) d\D°g\ if P 6 BV(Q) 



oo otherwise 



where in the case g G -BV^fi) the symbol -D a g denotes the Radon-Nikodym 
derivative of the absolute continuous part of the vector-valued Radon measure 
Dg with respect to Lebesgue measure and grgi^i is the Radon-Nikodym deriva- 
tive of the singular part D s g of Dg with respect to its total variation |-D S <?|. 
The function lo^ is the recession function: it is defined on fi, takes its values in 
the extended real numbers, and is given by 

ui^tx) _ ( K ltN \x\ if p= 1, 



<(*)= lim 



-oo if p > 1. 



for all ie!l. 

It follows that for p 



pm _j%|0/| if/eBV(0), 
bC Li(^)^UJ-| +oc if/ e 

and for p > 1 

qr - pm f^vl/II P if/e^(0), 
sc LHnrv)-\ +oo if/e L 1 (r>)\w rl - P (fi). 

Thus, the r^^Timit of the sequence of functionals 7?.^ has been established 
for all p > 1. By [5], example 1.21., it is clear that the first summand of the 
functionals is lower semicontinuous on L 1 (J7), and is thus -converging 
to the first summand of the functionals T v . Altogether, the following theorem 
has been proved. 

Theorem 7. For all 1 < p < oo the sequence of functionals (T p ) converges in 
the r^w^- sense to the limit functional T p . 

Corollary 8. For all 1 < p < oo the limit function f of the (sub)sequence of 
minimizers (/„) of !F P is the unique minimum of the limit functional T v over 
L 1 (f2). The minimum f also belongs to the space H W l p if p > 1 and 

Ll(n)nBV(n) ifp = i. 

Proof. The limit function / is a minimum of the limit functional T p over L 1 (f2) 
due to the properties of T~ -convergence [5]. It belongs to the space L 2 (Q) since 
S(f) is finite. Its mean is zero since it is the limit of the sequence (/„) from the 
closed subspace of L p (0), and it belongs to W 1 ^ if p > 1 and BV{9) if 

p = 1 by the third assertion of Proposition UJ □ 



3 Numerical Minimization of the Energy Func- 
tionals 

In this section we present some numerical schemes for minimiziation of the 
regularization functionals 7£„ in space dimensions one and two. We use a finite 
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element approach and approximate functions by a linear combination of finite 
elements, in particular piecewise constant and piecewise linear functions. Then, 
by tuning the kernel functions ip n we are able to recover standard finite difference 
schemes for total variation minimization on the one hand and on the other hand 
novel discrete schemes. A numerical comparison of the derived methods and 
applications to imaging will be studied in a forthcoming paper. The numerical 
schemes are derived from an approximation of the total variation functional in 
an infinite dimensional setting, which are then discretized. These schemes can 
serve as alternatives to existing numerical schemes (see e.g. [3 [TJ [TU1 H3 ITT] ) 
which are based on direct minimization of the Rudin-Osher-Fatemi functional 
and not on dual formulations, like the Chambolle's algorithm; see e.g. [4[ [2] . 



3.1 The One-dimensional Case 

We work on the domain 

Q:= (0,1). 

and consider minimization of TZ^ with a finite element method. 

1. The first two schemes are for piecewise constant finite elements: 

(a) We use the sequence of kernel functions (<£>„) defined by 

n 

'fin ■= -X[_i,i]. 

Let <xi,...,a n € R. Evaluating the one-dimensional piecewise con- 
stant function 



fn : = 22aiX[i=L,±] 
i=l 

with yields the standard Ty-seminorm of /„: 

n 

ni(f n ) = J2\ai-ai-i\ = \Df n \ 

i=2 

,(2)x 



(b) Using instead of (ip n ) the family of kernels (ip n ) defined by 

(21 n 



yields 

1 - 

2 



i=2 



We recall that ln(2) rj 0.7. 

2. Now, we consider a finite element method for piecewise linear splines. Let 
ao,...,a n € R and f n be the piecewise linear spline interpolating the 
nodes (|, a k ), k = 0...n, i.e., 

n 

fn := X! ai 9i 

i=0 
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where 



gi(x) := max I 1 — n 



x 

n 



.0 



Inserting /„ in 1Z\ and using the kernel functions ip n yields 



K-l(fn) = ^ 11 + t(ai- U ai, a l+1 ) 



n-1 



i=l 



i=l 



where 



|a ' +1 4 a ' 11 if sgn(ai_i - a,i) = sgn(ai - a i+1 ), 



We used the computer algebra program MAPLE for this evaluation. We 
provide a sketch of some parts of the computation for n>2. We want to 
evaluate 



10 jo 
which is equal to 

2 ^ L-i 



-<p n {x- y) dxdy. 



1(liB+i) IELo^)-EIU^)l dxdy 



F-y| 



Looking at the supports of the functions gi we realize that the latter double 
integral equals 



n 
2 



o Jo 



v+i 



ELo a i9i{x) - ELo a j9j(y) 



\x - y\ 



dxdy 



™ _1 ri rv+i 



k=2 



£"=n-2 <H9i{x) - E"=n-1 



dxdy+ 



\x-y\ 



dxdy 



We only treat the second double integral, the other two are evaluated 
analogously. Let 2 < k < n — 1 . 

Then 

rv+i 



k - y\ 



can be decomposed into the sum 



EiU-2 a%9i{x) - Ej=fe-i 



dxdy 



dxdy+ 
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k k 

r r 



EiU-i a t g l {x) - EjU-i aj9j(y) 



\x - 2/1 



x-y 



dxdy+ 



dxdy. 



whose summands we denote by Ik-, Ik and Ik,+ , respectively. The inte- 
grand of Ik is of the simple form 

| Oft - £Jfc_l| 

h = 



n 

By an application of Fubini's theorem 

h- = h-\,+- 

Thus, it suffices to treat the evaluation of Ik,- whose integrand J(x, y) is 
reshaped as follows: 

|(a fc -i - a fc )m/+(a fc -i - a fc _ 2 )nx+(a fc -2 + a fc - 2a fc _i)fc-a fc _ 2 + 2q fc _i-a fc | 

y - x 

From this representation it is already visible that the evaluation will be 
dependant from the sign of the differences 

A fc :=a fe -a fe _i and A fe _i = a fe _i - a fc _ 2 . 

We treat here one instance of the more complex case when 

sgn(A fe ) ^ sgn(A fe _i), 

namely the subcase where the middle coefficient a k -\ is the maximum 
of the three coefficients; the other subcase where ak-i is the minimum 
can be treated just the same. The less complex case where a k -\ lies 
between a k -2 and ak needs fewer case distinctions but apart from that 
can be treated analogously. Note that in the chosen subcase the second 
difference A^ = ak-2 — 2an + a k is negative. 

The numerator Jjy(x,y) of J(x,y) now reads 

J N (x,y) = \-A k ny + A k ^nx + A 2 k (k - 1)| 
and is positive if and only if 

A k ny- A 2 k (k- 1) 



x > fN(y) ■-- 



A fe _in 



Thus, to evaluate the inner integral of I k .- we have to determine the 
intersections of its integration domain (y — ^,-^-) with the intervals 
(/w(y),+oo) and (-oo, f N (y)), respectively, for all y e (^-, |), which is 
the integration domain of the outer integral. 

We get that 

, , , fc-1 k- 1 

f N (y) < < y 

n n 
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and 

y- l -<f N {y) <=► y< l{ k -% 

Let C := I (k - f|). Note that C < |. Then 

/■C* rfN{y) pC r^7T~ 

h- = / J{x,y)dxdy+\ / J{x,y) dxdy+ 

k k-1 

/ / J{x,y)dxdy. 

JC Jy-± 

Integrating J(x, y) with respect to x yields the primitive function 
K(x,y) := ((ra/- k + l)Al)ln(y - x) -mA k _i. 

Inserting the limits of the inner integral of the third summand yields 
K^y-^,yj = ((k-l-ny)Al) ln(n) - (ny - l)A fc _i, 

: " ((«y- fc + 1 ) A fc) ln (y-^) -(fc-i)A fc _i 



and 



L 3 (y):= / " J(s,y)da;= ((ny-A + l)Ag)ln(ny-fe + l)+(ni/-A;)A fc _i. 
A primitive function for L3 is 

(mi — k 4- 1 A 2 1 o 1 

M 3 (y) := ^ ^ ' k Hny-k+l)-—(ny - k + if A 2 k + -y(nk-2y)A k ^. 

Analogous computations for the first and second summand yield the func- 
tions 



L 2 (y) = (ny + l-k) (a 2 In (-^J 1 ) + A, 



and 

M 2 (y) = \ y(ny -2k + 2) U\ In ) + A, 



and 



L l{ y) = (-A 2 ) (ny-k + 1) (in (- A ^~* + *> ) - l) + A fe _ 



and 



. , A? / . „ „,., / A fe _i\ (ny- fc + l) 2 ln(m/- k + 1) 
Mi{y) = --± ly(-ny-2 + 2k)lnl—^)+^ L—AJ^ i 
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(ny — k + l) 2 + 2 ny (ny — 2k) 
2^ 



+ (ak - dfe-i) y, 



respectively. 
Now 



h, 



Mi 



k- 1 



Mi(C)+M 2 (C9-M2 



fc- 1 



M 3 



M 3 (C) 



(Ofc-l - Qfc) 2 + (Ofc-1 - afc-2) 2 



From this result and the results from above the final result follows easily. 
We further evaluated TZ^ for the Haar-functions using the kernel func- 
tions {(fin)- For j = k = the function is defined by 

for all ie!!. For k e N , 1 < j < 2 k we have 

( V¥ if T <= 

-V^ if a; €(1^,^), 
otherwise. 



Since /ig°^ is constant 



For we get 



ft 



I (4 0) ) - 



0. 



ft 



21n(2) ifn=l, 
2 ifn>l. 



Note that by symmetry for k > 1 and 1 < j < 2 



fti =fti (V fe 



fe-i 



(8) 



For fc > 1 the marginal functions h[ k ^ evaluate to 



ft 



' V2 r ((fc+2^r)M2)-(l-^)ln(2 fc -l)) if n = 1, 

^=((fc + 2)ln(2) -ln(n) + l) if 2 < n < 2 fe , 

V^n ln(2) - ^ ln(n) + - I) if 2 fe < n < 2 fc+1 , 

3V¥ ifn>2 fc+1 . 
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For k > 2 and j = 2, . . . , 2 k ~ 1 the inner functions hf- , evaluate to 



K 



i ("f ) 



/2fe 



f j ln(j) _ (j-l)ln(j-l) 
^ 2 fc 2 fc 

1 - In (2 k - j) + 

l-^)ln(2 fc -j + l) + ^) ifn=l, 



^(jln(i) + (i-l)ln(i-l) 
(fc + 2) In (2) - In (n) + 1) 



Ofc f (j-l)ln(j-l) 



) _ (j + l) ln(n) 
P 

(fcj+fe+2) ln(2) , 2+1 _ I 

2 k ' 2 k n , 



^((fc + l)ln(2)-ln(n) + l) 



if 2^ < „ < £ 
2 k_ J — — j ' 



ifC<n< 2 



j-i' 



if j 2 !! < n < 2fc+1 ) 



if n> 2 fc+1 . 



The evaluation of ft. ^ for j = 2 k 1 + 1 , . . . , 2 k is reduced to the evaluations 
directly above via ©. 



3.2 The Two-dimensional Case 

We now switch to the two-dimensional case and evaluate 71^ for a piecewise 
constant function defined on a subset of R 2 . In detail, let O be chosen as the 
open square (0, 1) x (0, 1) and / be defined on £1 via 

n 

for all ?; € f2 where := 7 fe x 7j with 7 fe := (^-, for all k, I € {1, . . . , n}. 
1. We choose the kernel functions 



<^n(u) := — X B ( ,i)(^) 

for all v in M 2 where 5 (0, ~) denotes the ball around the origin with radius 
— . The sequence (<p n ) satisfies all conditions stated in the introduction. 
We note that for two points (x, y) and (w,z) from O we have (x,y) — 
(w, z) € B (0, i) if and only if 



(w, z) € S u [ x, y, - ) := ((0, 1) x (0, 1)) n B I (x, y) 



We further define the intersection of the circle B((x,y), —) with the square 
h,l by 

Sk,i I x, y, - ) := 4,; H B I (x, y), 



18 



We have to evaluate 



K(f)= I I lf !f' V \ f ( {W '"h ((x,y) - (w,z)) d(w,z)d(x,y) 
JnJn \{x,y)-{w,z)\ 

n 2 f 1 f 1 ff E"j=i a i,jX/ M (a:,2/)-Efe,i=i a MX/ fe , i (w^) 
= — / / / / [7 m d ( w , z ) d V dx - 

7T Jo J0 JJSn(x,y i) \{X-W,y-Z)\ 

The occurring quadruple integral can be rewritten as follows: 

ttltlLli , ., ,,1? «»■'>«*■ 



i,.j = l fc,Z=l 



s fcii (x, y i)!^- w ,y- Z )\ 



For fixed 1 < i, j < n and fixed x,y G Iij the domain S^.; ^) of 
the inner double integral is empty if \i — fc| > 2 or |j — Z| > 2. Thus, it 
suffices to evaluate those summmands of the inner sum above that fulfill 
\i — k\ < 1 and \j — l\ < 1. However, in the case where i — k and j = I 
the integrand of the corresponding summand vanishes such that this case 
may be left out, too. 

Given a pair of indices (i, j) let 



{(k,l)€ {!,..., n} 2 : |i-fc|<l,|j-/|<l,(M)^(M)}- 



denote the set of pairs of indices for which the corresponding summands 
of the inner sum in the integral expression above do not vanish generally. 

Then the above quadruple integral equals 

Y] Y] \ai,j~ak,l\ [ [ (( 77 d(w,z)dydx. 

(ktx, , J^J^JJs^y^Kx ~w,y-z)\ 



We denote for all 1 < i, j < n and all (k,l) e Iij the quadruple integral 

jk,l 



expression on the right hand side above with J i '■ , i.e. 



J U ■= I I II 17 " r;d{w,z)dydx. 



,(x,y r)\{x-w,y- z)\ 



Again, let 1 < i,j < n be fixed. The set of pairs of indices with non- 
vanishing summands Tij may be partitioned into the sets 

2?j :={(M)e Iij : fc -and/ , ./! 

of pairs of indices marking squares diagonally adjacent to the square Iij 
and 

llj := {(k, I) G Iij : k = iorl=j} 

collecting pairs of indices that denote squares laterally adjacent to Iij. 
By simple transformations of the kind (x, w) \— ► (x ± — , w ± — ) etc. and 
applications of Fubini's theorem (as in the one-dimensional case) (or by 
geometric insight) we realize that the integrals J i ^ are equal for all (k, I) <E 
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ifj and the same holds true for all (k,l) € 2| •. Further, the respective 
values of the two evaluations are independent of i and j. 

Thus, it suffices to compute the values of Jl*^' J *~ and Jl'jf 1 for some 
fixed 2 < i t , < n, and the final result will be 

2 2 / / n n n—l n 

n n(f) = -J- 1 1 ~ a '- i .j- i i + 2 Z i aj j ~ °<+ij-ii ] J i.*j. 1 ' J * 

\\i=2 3=2 i=l j=2 



S 2 ~ a * j- 1 i + EE \ ai <o ~ j i J l*,l* 1 • 

K i=l j=2 i=2 j=l / / 

We begin with the more complex case of laterally adjacent squares and 
evaluate 

Jiff 1 = j j [[ d(w,z)dydx. 

h3 Ji^.Ji^.JJs iij _ 1 (x,v^)\(x -w,y-z)\ 

for some fixed 2 < i,j < n. 
We first point out that 



h:= 1 ^J^JL-^)\(*-™,y-*)\ d{w ' z)dydx 



such that 



^ - 1 rr 

d(w,z)dydx 

s tlj - 1 (x,v,i i )\(.x-'>»,V- z )\ 
J^]- 1 = 2Ji 

1 iJ 2 



This can be established by application of the transformations {x, w) i— » 

Let (x, y) be chosen from J^j with x > + We analyze the inner 
double integral 

// T7 " ^ d ( w , z ) 

JJs t , i - 1 (x,v,$)\{x-w,y-z)\ 

of Ji . We use the abbreviations 

2 

a := x — -, b := y — - and d := x — — . 

n n n 

for the distances of x and y to some nodes. Note that by our choice of x 
and y we have the inequalities a, b > 0, d < and, in particular, a > 

Let (w, z) e Sij-i (x,y, -). It follows that 

(x- W ) 2 < -L-(y-z) 2 . 



20 



and therefore, 



z > y and x - \ — - (y - z) 2 < w < x + \ — - (y - z) 2 . 

n | r V n 

Thus, z G (y — — , "^-) and we have to analyze the intersection of intervals 



l—,-)n{x-\l^-(y-z) 2 ^ + \l^-(y-z) 2 )- (io) 



(The index w in the symbol I w is just used as a symbol to indicate that 
we are dealing with the integration domain of the variable w but does not 
stand for the values of w. The same applies to I z etc. below.) We first 
point out that 



i-l I 1 / 1 

72 



< x - \l ~r - (y - z ) 2 <== ^ z < y - \ -? - a2 - ( n ) 



n V n VTi 

Thus the result of ifTO)) is dependent from the intersection 



1 =[y-^ J —^) n (-^y-\l^~ a2 • (12) 



nr 



While it is clear that y — J X — a 2 > y — — we have 



y-\\- 2 -a 2 < 3 — <=> y<^ + J^-a 2 . I IS) 
" n J n n V ^ 



with i=i < + JX -a 2 <i- 



We first consider the case that y G I y := ( 2 — , ^— + J — a 2 j . Then 



h = (y-hy-^-a 2) 

We look at the subcase z G I z . It follows by (JTTJ) that now the lower 
bound of I w is x - ^ - {y - z) 2 . 

Considering its upper bound we have to find the minimum of 
x + \ I X — {y — z) 2 and — . Similarly to (|TT|) we get that 



- > x + J-s -{y-zf ^ z<y- J— - d 2 . (14) 



Since a > d it is true that y — J X — d 2 G I z implying that the treated 



subcase has two more subsubcases: z G I z ,i '■= ( y — z:, y — \ — — d 2 



and z G I z ,2 ■= [y~ y ~d 2 ,y- 

By ([IH and (fTU|) the upper bound of !„, is x + J 4^ — (y — z) 2 if z G I 2 ,i 
and equals — if z G I z> %. Thus, the treated subcase gives rise to the 
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following two quadruple integrals 



dwdzdydx 



v /j ? _ fo _ z) 2 \(x-w,y-z)\ 
1 



=1 Jy-J^TZ-p J x _/^_ iy _ z y \{x-w,y-z 

which we denote by _ftTi and K2, respectively. 



dwdzdydx 



We turn to the subcase that z £ I z , i.e. z G ( y — i /-A — a 2 , ■ 2 — ) . By 



(fTTj) the lower bound of I w is in this subcase . Since y — J — a 2 > 



V - y ^2 - d 2 by (H3J) in this subcase the upper bound of I w is ^. This 
subcase yields the quadruple integral 



K 3 -= I I / , / J7 rj dwdzdydx. 

Jy-^-adi^ \{x-w,y- z)\ 



We still need to analyze the case y £ I y , that is, y £ {^~7T + \j ~ ° 2 ' n) 
In this case by (fT5|) and (fT^)) the interval I z equals (y — ^-).By JT 



and it is clear that the lower bound of I w is x — y 4y — (y — z) 2 . The 
determination of the upper bound is a little more intricate including two 
subcases concerning the choice of the domain of y one of which generating 
two subsubcases concerning the domain of z. However, its computation 
processes similarly enough to the computations in the first case that we 
skip it here and just state the resulting quadruple integrals which we name 
K4,K§ and Kg, respectively. 



dwdzdydx. 



J y -± 7 a; - v /Jj-( 2/ - z )2 \[x-w,y-z 

dwdzdydx 



i__i + j_7i_i +v /^^ y i/ _^rz^i x _ v ^rz7— j? |(x- z)| 



'i^l+i Ji^+JjTZ^Jy-i- J x -^-^ z )2 \(x-w,y-z)\ 
Altogether, 



dwdzdydx. 



The evaluation of the six quadruple integrals if^ involves the transfor- 
mation of the respective inner double integrals to polar coordinates. In 
order to simplify this procedure we first translate the integration domain 
of the respective inner double integral to the rectangle (— — , — ) x (0, -). 
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In all six cases given a point (x, y) from the domain of the respective 
outer double integral this is done by application of the transformation 
(w, z) i ► (x — w,y — z). Let Li be the result of this application to Ki. 
Then 



Li = / / / j- -7 dwdzdydx, 

7i_i +5 L7i_i JjJ—Z2j-^—2 \{w,z)\ 

ri r^+^- a2 /V^-" 2 /V^ 2 1 

L 2 = / / / ,7 ry dwdzdydx, 

Jl^L + ^Jl^. J d \{W,Z)\ 



L 3 = / / / / dwdzdydx, 

Ji=l + iji=l J b J d \(w,z)\ 



ft f^r+\J-^- a ' fn f\j-^" z ' 1 
Li = / / / 77 ry dwdzdydx, 



/ / / / 77 77 dwdzdydx, 

J i^ + ^J i=± + jj--; J b J d \(w,zy 

i 3 1 / rT T 

r r n 



dwdzdydx. 



Let for 1 < i < 6 the function Fi be the evaluation function of the inner 
double integral of Li defined on the domain of the outer double integral 
of Li. 

In Li,L 4 and L§ the integration domain of the inner double integral is 
a segment of the circle B(0, i) that results from the intersection of that 
circle with a parallel to the x-axis. A straight-forward transformation to 
polar coordinates (r, <fi) yields for example for F\ : 

- V " r I +7T 

*i(s,J/) = /, / /at— *x ' d<j>dr 



dr 




= n r — 2 r arcsin 



^— — d? arcoth 



7T arCSm {\[^~^ n ) 

n n 



- rf 2 ( In I - + d ] - In ( - - d 



n 2 \ \n I \n 
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and F4, and Fq are treated analogously. 

The integration domain of F3 is a rectangle with edges parallel to the 
axes stretching across both quadrants of the upper half plane. In order to 
transform this domain to polar coordinates we split it along the y-axis in 
two axis-parallel rectangles that reside in the second and first quadrant, 
respectively, 



F 3 (x,y) = 



1 



\(w,z)\ 



dwdz 



y5" 



1 



IO, z )l 



dwdz, 



and call the resulting double integrals A(x,y) and B(x,y). 

We turn to the computation of A(x, y). The transformation of a rectangle 
domain located in the second quadrant to polar coordinates depends on 
whether its bottom left vertex or its top right vertex is more distant from 
the origin. In the case of A(x,y) this conditions reads 



\(d,b)\ < 



0, 



1 



(15) 



In the case where (1151) holds true the transformation to polar coordinates 
yields 

p\(d,b)\ p— arcsin(£)+7r p | (o, 4j -a 2 ) | /-arccos(|) 

A(x, y) = / / d(j)dr + / / dcj)dr+ 

Jb Jf J\(d,b)\ 



(O,^-— a 2 )| J — arcsin ^ 1 r J +tt 

in the opposite case A{x, y) equals 

|(0,^--a 2 )| p- arcsin(£)+7r p\(d-b)\ p- arcsin(£ )+7r 

/ d(f>dr+ / / / 2 \ d(pdr+ 

(^d,^J~-X 7 — a 2 ^j | />arccos(^) 



J — arcsin I * ^ +7r 



In both cases the occurring three double integrals can be evaluated sim- 
ilarly like Fi(x, y) above. Summing together the respective three results 
yields in both cases the same result: a sum consisting of summands that 
are of one of the following three types: binary products where one factor 
is a logarithmic expression, binary products where one factor is an arcsin- 
or arccos-expression or binary products of a square root and ir. By use of 
the appropriate transformation rules for arcus-expressions the two latter 
groups of binary products cancel each other out. Therefore, in both cases 
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A(x,y) equals 

^&(ln(Vd 2 + & 2 +d) - ln(y/cP + lP -d)) + 

d(^ In i^s/d 2 + b 2 + bj - In i^s/d 2 + b 2 -bj + 




1 - a 2 + d 2 - d | - In ( W 4 - a 2 + d 2 + d 



n 2 \ \ n 2 



The double integral B(x, y) is evaluated in a completely analogous fashion. 

The integration domains of F 2 and F 5 have a similar geometric structure. 
We give a short overview of the evaluation of F^- As in the case of F$ 
we split the integration domain along the j/-axis in order to have less case 
distinction when transforming to polar coordinates: 



f\J^- d2 f° 1 f\J^- d ' 2 f\[^- z2 1 

F2(x,y) := / / j- ry dwdz+ / / dydx. 

JJjTZ^ J d \{w,z)\ Jjj—; Jv \{w,z)\ 



The first double integral has an axis-parallel rectangle domain located in 
the second quadrant and is treated like A(x, y). The domain of the second 
double integral is the intersection of two shapes: an axis-parallel rectangle 
domain located in the first quadrant whose lower right vertice B lies on 
the circle B(0, -), and the circle B(0, — ) itself. This means that by the 
intersection the right edge of the rectangle and parts of its top edge are 
exchanged with a circular arc around zero with radius ^. The double 
integral is translated to polar coordinates as follows: 



r\J-^- d2 rt rt 



d(f>dr. 



For the evaluation of F$ proceed as for F2 . The only difference to F2 lies 
in the fact that by the intersection with the circle B(0, i) also parts of 
the bottom line of the corresponding underlying rectangle are removed. 
In symbols this is reflected by exchanging every occurrence of the term 



— a 2 with b in the double integral directly above. 

The integration of the functions Fj with respect to y can be executed in 
all cases by standard means. Note that the double integrals related to 
Fi, F2 and F4 do not depend on y such that those integrations are mere 
multiplications of the respective functions with the difference between the 
limits of the respective integrals. After all the resulting functions of the 
variable x have been summed up the fourth integration can be carried out 
yielding 

6 



2 



3n 3 * J 5 3„3' 
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This solves the case of laterally adjacent squares 

The case of diagonally adjacent squares is much simpler. For 2 < i, j < n 
we have to compute 

T i-i,j-i _ f n f n f f 1 



J = / / / 77 tj d(w, z) dydx. 

where 

5 — {*>y>k) = ((^T'^r) x (^^)) n K (x '^ 

An easy computation shows that the latter set is empty if and only if 
| (a, 6)| > i. Therefore, for the computation of the inner double integral 
of J?~~ we may restrict ourselves to points {x,y) € Ii.j that satisfy 
| (a, 6)| < —. By similar reasoning as in the lateral case we infer that 
J*~ lj_1 equals 



s /• 2 ^ i +v / ^- a2 



dwdzdydx. 



By applying the transformation (w, z) i— > (x — w,y — z) this quadruple 
integral transforms to 



^ + ^-« 2 fy/i*- a * [\Jt?- z2 1 

dwdzdydx, 



lb J a \{w,z)\ 

and by changing the inner double integral to polar coordinates we get 



r^<^ Jr \J ~^~ a2 fn /-arccos^) 

J\~ 3 = I I II — 77 dwdzdydx 



-i /j-i 

I 

6n 3 



|(a,6)| Arcsm(i) K^ 2 )! 



Altogether, by (O the final result is 

x n n n— 1 n 



i=2 3=2 i=l 3=2 



.inn n n 

yi=l j=2 4 =2 3=1 

2. With / as above we also evaluated TZ^f) with the kernel functions 

o 

n 

<Pn{v) := — i)x(_ii)(«) 



2G 



for all v in E 2 . Note that these kernel functions are not radial such that 
this case is not treated in the theory section. The evaluation proceeds 
similar to the above one. As a result we get 

1 ^f^2 1 / 71 U n—l n 

T^n(f) = 3 H " + H ^ ~ 

\i=2 j=2 i=l j=2 

1 3 In ( v / 2 + 1) - 3 In (v^- 1) - 2 (a/2 - 1) 
12 n X 

n ii n n 

X! H i a '-j ~ "• .• 1 i ai J ~ j i 

v i=l j=2 i=2 j = l 
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